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ABSTRACT 



I 



The Coppens-Sanders theoary for the one-dimensional, nonlinear, 
acoustic wave equation with dissipative term describing the viscous 
and thermal energy losses encoiantered in a rigid-walled, closed tube 
of large length-to-diameter ratio was applied to finite-aiiplitude 
standing waves by the use of the Fast Fourier Transform. Conputer 
programs were written to determine the amplitudes and phases of the 
first 255 harmonics. Curves of harmonic distortion as a function of 
the strength p>arameter were foiand to be in excellent agreement with 
available exp>erimental data, to agree with the Coppjens-Sanders 
pierturbation analysis, and to extend the theoretically describable 
regime closer to strengths leading to the formation of the shock front. 
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r 

phase speed in the tube associated with the n^ harmonic 
1/2 

^9p/8p) ' _ for dissipationless propagation (linearized theory) 



co/c 

length of the tube 

Lagrangian coordinate measured in the x direction frcan the 
driven end 

particle displacorent in the x direction 
particle speed in the x direction, (9C/8t) 
acoustic pressure 

acoustic pressxjre for the n^ harmonic 
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2 

P /Mp c 
n o o 

phase angle associated with the tonporal part of the acoustic 
pressure 



dissipation parameter 



th 



dissipation parameter for the n frecjuency conponent 

D 'Alambertian in Lagrangian coordinates (one spatial dimension) 

th 

dissipation operator for the n^ frequency conponent 



(6 /co )sV9^3t - 6 3^/9>c^ 
n n ^ X n ' 

a Mach number (See Eq. 2.3) 

2 2 

parameter of nonlinearity p (3 p/3p ) 



b =1+1/2 (B/A) 

= amplitude operator (See Eq. 2.13) 

9^ = phase operator (See Eq. 2.14) 

2 

= amplitude associated with p^ (See Eq. 2.9) 
= phase angle associated with B^ 

Hg = coefficient of bulk viscosity 
n = coefficient of shear viscosity 
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I. INTRODUCTION 



Theoretical investigations of finite-amplitude waves have been 
characterized by various approximation and perturbation techniques, 
but more significantly by the postulation of various models of the 
absorptive mechanisriB [1,2,3,4] . Coppens and Sanders [5] have in- 
vestigated finite-amplitude standing waves for rigid-walled tubes 
of large length-to-diameter ratio by a perturbation method and have 
proposed that the doidnant energy losses are produced by the viscous 
effects of the fluid at the walls and the thermal losses due to the 
wall material. They have shown that an approximate wave equation for 



the process is 



u 



J_ \ ^ I 



\ 2 . 



c)x 



( 1 . 1 ) 



n-\ 

where LL 



CO 
ST ^ 



th 



and the subscript n represents the n^‘ harmonic 
of the waveform. The n^ frequency carponent of<2) is 









8 ^ 



( 1 . 2 ) 



th 

where 6^ is the dissipation parameter [3,5] for the n fre(guency 
conponent defined by = '^1^'Vn* process has two characteristic 

parameters, one related to the strength of the standing wave and the 
other related to its frequency. 

The purpose of this thesis is extend the Coppens-Sanders approacdr 
by a "Fast Fourier" analysis technique. By an iterative method, the 
amplitudes and phases of the various harmonics are obtained for dif- 
ferent values of the strength and frequen(oy parameters and are conpared 
with available experimental and theoretical results. 
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II . THEORY 



In this section a Fourier-synthesis approach is presented for 
the one-dimensional, nonlinear, dissipative, acoustic wave equation 
based on the theory of Coppens and Sanders as applied to finite- 
amplitude standing waves in a rigid-walled, closed tube of large 
length-to-diameter ratio. 

We begin with Eq. 1.1 which describes the process in question. 
Notice that 



oo 









( 2 . 1 ) 



so that Eq. 1.1 may be written as 



OO 

n-\ 













2- 






( 2 . 2 ) 



Assume that p may be written as a Fourier series 
oo 

COSX\Jk[L- y) (ncjt+ I (2.3) 

n-\ 

2 

where P = Mp c A . 
n o o n 

2 

Substitution forO.. and0r into the left-hand side of Eq. 2.2 

L Ln 

results in 

OO 



n=\ 



S ^ 






lPn‘ (2.4) 



Upon substitution for 6 and w , Exp. 2.4 becomes 
^ n n 

L u_ S, 

^ n . o "> 



S \ o 



V\ 



(2.5) 



a=\ 



14 



Notice that ^ — _ 



ni^lVn 

/ oX 



noj 



■(Pn 



SO that Exp. 2.5 may be written as 
Co 



VinSk 



r\= \ 



2 



S, 






■f 



co 



S, S 






n 



( 2 . 6 ) 



Define w =co + Aoj where Aw << oj , and for the fundamental = 

^ ^ C, 

viiere the approximation is valid for weak absorption. The approximation 



may be made that 



CJ 



2 _ 2 
o 






e. 



+ 






2Aoj 



C, 



CJ, 



(2.7) 



Define (c^/c^) = 1 - 6 ^^ and neglect second-order terms; Exp. 2.6 may 



then be written as 

^ .2-r . 



-SV (n5^) 



n = \ 



.n' 



ZAcj 

t^rS, 



i 









n 



( 2 . 8 ) 



The right-hand side of Eq. 2.2 contains the quantity p . This 
involves multiplication of two infinite series. Terms of the form 
cos p(L-x) sin (qwt + <|) ) for p 7 ^ q may be neglected [5 ] , so that p"^ 



may be approxinnated by 

B.cosni^ 



L-X 



SWl 



ncij 






(2.9) 



'TX- \ 



Then we have tm- written as 



00 



s. 






Ijn^ 



ZAcj _j ^ 

S, 



AnCosn.fejL-X ^ s\njnajt-^j(2 • lo) 



00 



_Mb 
- 2 . 




n 



cosn 



9, 



L-x 



s\n ncj 



A-t-r 



n-\ 
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r|=i 



If attention is limited to the rigid end at x = L, Eq. 2.10 may be 

considerably simplified to 
OO 

2A0J 






Ln 



■/2 






Si ^ r - 



This becones 

OO 



<oo 

— ^ ^ ^ Sri s'\ n (v^cjt * 



n=\ 



V An‘Sm(noJt 

OO 



( 2 . 12 ) 






v^ere 






and 



L\n"^ ' + ^r 



“rf- ' 






n 



'A 



/ J 



(2.13) 



^n-\ 



-I 

an 



n 



'/2 



Vz 



ZAcJ 



(2.14) 



n " ”r^i 

The parameters relating to the strength and frequency response of the 

standing wave are Mb/6j^ and , respectively. 

r 1 

An iterative solution to Eq. 1.1 may be obtained by successive 
application of Eqs. 2.12, 2.13, and 2.14 to an initial set of A^'s 
and The necessary procedure is shown as a flow diagram in Fig. 

1. A discussion of the Fast Fourier Transform and the reasons for its 
use herein are found in Appendix A. Notice that any D.C. terms which 
appear in the procedure have to be set equal to zero since the theory 
excludes these terms. 
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A soiarce of cxanputational difficulty arises vAien the frequency 

2A(jl) 

parameter ^ beccmes unity. For this case Eq. 2.13 reduces 

“r 1 

to =\/n/2^ a result which causes the ccmputer solution to become 

unstable because of the strong excitation of the higher harmonics. 

The reason for this difficulty may be seen as follows: the resonance 

frequencies for the classical nondissipative case are defined by 

0 ) = nu , where a> = c . For linearized propagation in a duct, 

resonance frequencies for exciting each harmonic would be given by 

' = no>' = n^ c , where c = c (1-6 /2) . As n we have the 
n Tj n no n 

« 

limits c -> c and a>' = w . For the case of large n then, since 
n o c ^ ' 

u) = + Aw, the value of the response parameter corresponding to the 

, 2 Aoj 

predicted resonances in the linearized theory is given by ■ = 1; 

r 1 

thus, the higher harmonics are optimally excited and problems will 
arise in the ccsiputer program, vhich admits only the first 255 harmonics. 

Notice that the bulk losses have been neglected up to this point. 

If these terms are retainedO^ may be rewritten [5] as 



(D, - 

Ln — 



s 



D_ 



(jJ 









z 



-sJ 






(2.15) 



In most experimental configurations the value of e is usually so small 
that it may be neglected when compared to Injection of an 

unrealistically large value of e does alleviate, to sane extent, the 



instabilities in the neighborhood of 
and 2.14 becane 



2Aw 

w 6, 

r 1 



= 1 



Equations 2.13 






Lin''" 



4* 



2Acd 



2 






r / 



2 



n 



77 - + € n 

Va / 



( 2 . 16 ) 
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and 





-€n 




n 



Vz 






2A(>l) 

^v-S, 



(2.17) 



f =^(poCcf)'(^S+ 3-’?) - 

and would have the asymptotic value of 1/e n as n The 

effectiveness of this modification will be discussed in more detail 
in Section III. 
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Flow Diagram for Solving Eq. 1.1 
Figure 1 
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III. RESULTS, DISCUSSION, AISD CONCLUSIONS 



The computer program written to solve Eq. 1.1 using Eqs. 2.12, 

2.13, and 2.14 gave results consistent with the theoretical pre- 
dictions of Coppens and Sanders [5] and Ruff [10] . An additional 
perturbation solution that was performed to insure that the program 
was operating properly may be found in Appendix B. An initially 
sinusoidal pressure waveform when used with a correction factor of 
unity gave rapid convergence to a distorted waveform except that in 

the neighborhood of response parameter — the solution 

•"r^l 

became unstable. This region of instability gradually increased as 
the strength parameter was increased. For these cases the modifications 
defined by Eqs. 2.16 and 2.17 were used with e = 0.0 (Hand allowed a few 
more values of to be ccmputed, but was not successful in gener- 

ating stable results for the higher strengths. The excellent agree- 
ment with Coppens and Sanders is illustrated in Fig. 2 for the strength 
curves. Figures 3, 4, and 5 illustrate both the area vAiere the solution 
deteriorated and also the frequency characteristics of the various 
harmonics. 

Since the only difficulties encoiantered were in the vicinity of a 
response parameter of lanity, it seems plausible that the instability is 
a result of the behavior of the higher harmonics which beocme in^rtant 
in this region. Recall that a response parameter of unity corresponds 
to the optimum excitation of the higher harmonics and that of necessity 
a truncated series has been used. Then the higher harmonics that are 
retained are affected by the absence of the harmonics that have been 
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discarded and react by graving in amplitude to absorb the energy tliat 
would have been dissipated by those terms that have been emitted. The 
next refinement of the theoretical development of this research would 
be to cause the retained harmonics to behave as if the neglected higher 
harmonics were present. 

Since a solution using iterative processes is easily adapted to a 
cemputer simulation it is suggested that this method be continued. The 
use of a time-sharing system is recommended since it allows close control 
of the program and saves a great amount of time. 

It has been shown that a Fourier-analysis approach to the Coppens- 
Sanders theory for the one-dimensional, nonlinear acoustic wave equation 
is readily adapted to cerrputer simulation and extends the region of 
validity to the pre-shock r^ime. The series truncation appears to be 
the only c±>stacle to obtaining a solution valid for the entire spectrum 
of strength and response p>arameters. 
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Strength Curves 
Figure 2 
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Q Withcait bulk absorption 
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Q Without bulk absorption 




u £. 

r 1 

Response Curves, Mb/6j^ = 0.4 
Figixre 4 
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CO 6-, 

r 1 

RespcDnse Curves, Mb/6j^ = 0.6 



Figure 5 




Response CurvM/ Mb/6j^ = 0.8 
Figure 6 



26 



APPENDIX A 



The Fast Foiarier Transform 

The Fast Fourier Transform [6, 7, 8, 9] provides an effective means 
for estimating the Fourier coefficients of a set of evenly-spaced 
discrete data. The Fast Fourier coefficients differ only slightly 
fran the conventional Fourier coefficients. Construction of wave- 
forms is of prime interest for the higher strength parameters, where 
the Fast Fourier Transform has the distinct advantage of eliminating 
the Gaussian overshoot characteristic of the conventional Fourier 
Transform. This property is illustrated in Figs. 7, 8, and 9. The 
waveform of Fig. 7 was sampled to provide 512 evenly-spaced discrete 
data. This was the input to the procedures available on the IBM-360 
Cdiputer that perform Fourier and Fast Fourier analyses. The coef- 
ficients frcm the Fourier analysis were used to construct the waveform 
of Fig. 8, while those of the Fast Fourier analysis were used for Fig. 
9. Fran the standpoint of computer utilization, the Fast Fourier 
Transform has the advantages of requiring less space and less time 
and was used for the analyses in this research. 
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Waveform for Analysis 
Figure 7 
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Reconstxucted Waveform Using Fourier Transform 

Figure 8 
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Reconstructed Waveform Using Fast Fourier Transform 

Figure 9 
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APPENDIX B 



Perturbation Solution of Eq. 1.1 

Assume as a first approximation that the solution to Eq. 1.1 
consists of only the fundamental which may be written as 
Z, 



- PoCq\\A COS . 



(B.l) 



Then 






PqCqM 

2 



cos 2 i^(L-x) -COSZcdt 

-COSzJ^(L-x) COSZcdt 



(B.2) 



so that the right-hand side of Eq. 2.2 is, at x = L, 



_b_ ZM)- hp„C ^Mj^^-|+COS2ajt] . 

ox 

The squaring of the sin((jjt) gives a D.C. term and cos(2ojt) . 

The D.C. terms may be excluded [5] so that the second harmonic required 
on the left-hand side of Eq. 2.2 may be assumed to be of the form 

p2^ COSZ^(L-x)s^'n(Zcjt . (B.4) 

/ 2 . 

Then - z') -6^ 2. may be written as 

^ Co y 2- ^ 

+ zVzT' I C0S2^ (z- P)coZ2t^t+<)5>2')®- 

2 

For Au) = 0, JiL — I , so that at x = L, Exp. B.5 may be written as 



-APjiifS, l-h/vz’ ]sin(2o)t-‘-<i>J-t4/T')cos(2wt+4) 



(B.6) 
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v^ich when simplified is 



.7G8s'\xn(zcot + c^)2 + 247.b° 



z 



(B.7) 



Equating Exps. B.3 and B.7 

§1 .7(oS S\n b°) 

2 2 2r 1 

-bp^c^ W\Sk L-l+eosZwt 

Mb 

which may be solved for P 2 and (j> 2 . For -^ = 0.1, 

.0525poCt^M and dz - . 89rad\an s . 

This agrees with the predicted values of Coppens and Sanders [5] 



and Ruff [10] and the results of the computer program written to solve 



Eq. 1.1. 



32 



COMPUTER PROGRAM 



This section contains the canputer program losed to solve Eq. 1.1 
and a list of the symbology pertinent to the program. 



AAMP(N) 
FEE(N) 
BAMP(N) 
GAMMA (N) 
H(N) 

THETA (N) 
P(N) 

R(N) 

ARG 

RESPM 

STRPM 

CORR 



B 



H 



e 



n-l 

^n-1 

n-l 

^n-1 

n-l 

n-l 



amplitudes of p/Mp^c^ for increment N 

2 2 

amplitudes of (p/Mp^c^ ) for increment N 
o)t/511 



2Aoj/ 

Mb/6j^ 

Correction factor 



Canputer Program Syn±ology 
Table I 
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Cotputer Program for Solvixig Eq. 1.1 



niMc^NKlTTN A^wp(2'=7),'^Ff'{2'^'’),R^'«0f2'^'') »r,* MMA ( 2R7 ), 
1H( 2K7) ,THFTA (2F7) , TMV(f4) »S (64) ,C(600 ) , « (600 ) 

PT=3. 141^0^ 

4pr, = 2*«’T /Fll . 



TNITIfl T7C PE:CDM,STppM,COPR,AAMDf N) tPCCfM), 



OPAO (6,2) OP SPM,FTRPM,cnpo 
2 FPPMt T( 3F6.3) 

OP.' 0(5,^) (MMPH) ,FCP (I ) ,T=1 ,2^^) 

4 COPMA T( 3(4X,2F1C,7) ) 



thip crcTIPM C^L('ULATFP />MCL(T|inF oopp^^TFR*^ H(P) AMT 

Ti^c PHAOp PTFOA'^PR" T^-FTi(M), 



no 1 1=2,266 

J = T-1 

XM=J 

06TIP=-( l.C/S00T(XM) ) 

ni)MMYl=-PA'^T n- 1 . 0 <-RFFO^^ 

CQUA'^inN 2.13. 

H( T ) = l. C/SnoO(0(|MMYl *n'IMWYH-RAT Tp«P AT ID 
EOUAtipa' 2.14. 

TMpTf ( I )=AT^ N2( PATin,CilMMYl ) 

1 C^N'TINliF 

O^ 1000 MPASC=1 ,27 



THIS ocpTjDM calculates TH^ AMplITUOF'P BAMP(N) ANO '^HF 
r>MASFP PAMMAfK) CCRRE^POMni MG TP TH= FQijapp op tmc 

popccyRO. 



11 1 = 1,612 
o( I )= 0 .C 
np 12 J=l ,257 

o( T ) = o( I )+AA MO( J)*SIN( ARG*(.J-1 )«(I-l) + EEo(J)) 

12 CHM'^INIIF 

R ( I )=( o( T )*D(I) ) 

11 COM'^TNUF 

ts TMPljT Q(I) CPMTAIKS TH= AMPLITUOF^, fqp INCREMEMT 
"T,'» CroRF SPOMPI MG Tp Tuc SOUARF PF THF 

CALL RM/RV ( R,P ,TMV,‘^ , TEFRO ) 

AS OUTOUT, P(I) CCNTAIMS Tmc eaRT FOURIER COEFCICIEMTS 
OF THE SQUARE OF TMp pqesSUP=: (A 0 )/ 2 , BC= 0 , Ai, Rl,.. 

(AN)/ 2 , BM= 0 . 

TWF am«c top TWP COEFEICTENTS OF THE CGSTNR TcpwR AMO 
THE BNtc aPF the COEFFICIENTS OF THE RIME TERMS. 

DO 13 1 = 3 , 511,2 

WPItf each CCMPOMEMT OE Tpr cqdapc qf the opF^cijoF AS 
RiMO Cl M( M( i MCRFMEMT) + gamma (N)) RY !|F IMG THE fact 

FPURTFR COFFcir Ic^Jp <5 , 

J=( T + 1 ) /2 

RAMP( J)=<0RT(R{I + 1)*P(H-1)+R(I)«R{ I) ) 
gamma ( J)=A taN2( P ( I ) »R (T + 1) ) 

13 CONTI nUF 
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SFT O.C. TPRVS TC 7 FRG. 

RAMD ( l )= C.O 
RAMDf 2'^7)=0.C 

GAMMA { 1 ) = 0.0 

gamma ( ? c 7 ) = 0.0 



THIf? SFCtIOK' adjusts the LFFT-HANT side of EOUA'^ION 2.12 

TOUAPnS THE right-hand SIDE. 



nn 22 1=3,256 

AAmp( T )=AAmp(X)+CORR*(BAMP{I)*STRDm«H( I )* 0 . 5-AAMP( I ) ) 
FCFf I )=FpF(i ) + (GAMMA(I )-pfh( I)-'^hcta(I ) )«C0RR 
22 COMTTNU'^ 

SPT D.r . TERMS TG 7ERC. 

AAMRI 1 )=C.n 
PFi=( 1 )=C.C 
AAMO( 257)=C.0 
FRC( 257)=C.O 



PPIN"^ AMD PUNCH CUTPU"^. 



34 



35 

32 



33 

31 

lOOD 



45 

43 



WRTtf (6,34) NPASS 
FnpMA.T(« pAS^ MUMRFR *,T5) 
nn 32 1=1,10 
J = I-1 



? IT? ( 6,35) J,AAMP(I ) ,Fi=F( I) ,B? MO( I) ,gamma( I ) 

IP MAT (I5,4F15.7) 

■>mtiNUF 

^ 31 1=11,257,10 
= 1-1 

>TTF (6,33) J,AAMP(I ) ,FFF( I) 
tPMATI !5,2F15.7) 

■' MtinijF 

•’NTT N(ir 

T'^F (6,45) RE SPM,STRPm,COPR 

>ma T( lOX, 'RE5PM= • ,P6.3 ,5X, •STRPM=* , F 5 . 3 , 5X , • C0RR= ' 

I 3 ) 



.3) 

MCH 43 ,(I ,AAMP( T) ,FF«=(I ) ,1=1 
5MAT(3( I4,2F10.7) ) 



257 ) 



END 



f 
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